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ABSTRACT 

We present a detailed analysis of week-long simultaneous observations of the blazar Mrk 421 at 2-60 keV 
X-rays {RossiXTE) and TeV 7-rays (Whipple and HEGRA) in 2001. Accompanying optical monitoring was 
performed with the Mt. Hopkins 48" telescope. The unprecedented quality of this dataset enables us to establish 
the existence of the correlation between the TeV and X-ray luminosities, and also to start unveiling some of its 
characteristics, in particular its energy dependence, and time variability. The source shows strong variations in 
both X-ray and 7-ray bands, which are highly correlated. No evidence of a X-ray /7-ray interband lag r is found 
on the full week dataset, with r < 3ks. A detailed analysis of the March 19 flare, however, reveals that data 
are not consistent with the peak of the outburst in the 2^keV X-ray and TeV band being simultaneous. We 
estimate a 2. 1 ± 0.7 ks TeV lag. The amplitudes of the X-ray and 7-ray variations are also highly correlated, and 
the TeV luminosity increases more than linearly with respect to the X-ray one. The high degree of correlation 
lends further support to the standard model in which a unique electrons population produces the X-rays by 
synchrotron radiation and the 7-ray component by inverse Compton scattering. However, the finding that for 
the individual best observed flares the 7-ray flux scales approximately quadratically with respect to the X-ray 
flux, poses a serious challenge to emission models for TeV blazars, as it requires rather special conditions 
and/or fine tuning of the temporal evolution of the physical parameters of the emission region. We briefly 
discuss the astrophysical consequences of these new findings in the context of the competing models for the jet 
emission in blazars. 

Subject headings: galaxies: active — galaxies: jets — BL Lacertae objects: individual (Mrk 421) — gamma- 
rays: observations — X-rays: individual (Mrk 42 1 ) — radiation mechanisms: non-thermal 
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1. INTRODUCTION 

Mrk 42 1 is the brightest BL Lac object in the X-ray and UV 
sky and the first extragalactic source detected at TeV energies 
(Punch et al. 1992). Like most blazars, its spectral energy 
distribution shows two smooth broad band components (e.g. 
Sambruna, Maraschi, & Urry 1996; Ulrich, Maraschi, & Urry 
1997; Fossati et al. 1998). The first one extends from radio to 
X-rays with a peak in the soft to medium X-ray range; the sec- 
ond one extends up to the GeV to TeV energies, with a peak 
presumed to be around 100 GeV. The emission up to X-rays is 
thought to be due to synchrotron radiation from high-energy 
electrons, while the origin of the luminous 7-ray radiation is 
more uncertain. Possibilities include inverse Compton scat- 
tering of synchrotron (synchro self-Compton, SSC) or am- 
bient photons (external Compton, EC) off a single electron 
population thus accounting for the spectral "similarity" of the 
two components (e.g. Macomb et al. 1995; Mastichiadis & 
Kirk 1997; Tavecchio et al. 1998; Maraschi et al. 1992; Sikora 
et al. 1994; Dermer et al. 1992). Alternative "hadronic" mod- 
els produce 7-ray from protons, either directly (proton syn- 
chrotron) or indirectly (e.g. synchrotron from a second elec- 
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tron population produced by a cascade induced by the inter- 
action of high-energy protons with ambient photons) (Miicke 
et al. 2003; Bottcher & Reimer 2004). The synchrotron pro- 
ton scenario may be more favorable for objects like Mrk 421 
(Miicke et al. 2003), because of the lower density of the dif- 
fuse photon fields necessary for processes like pion photo- 
production to be effective. Moreover, it is generally true that 
"hadronic" models need a higher level of tuning in order to re- 
produce the observed highly correlated X-ray/7-ray variabil- 
ity. Hence, in this paper we have not addressed this class of 
models, and instead focused our limited modeling effort on 
the pure SSC model. 

All of the above models of the 7-ray emission from blazars 
have all had some degree of success in reproducing both 
single-epoch spectral energy distributions and their relative 
epoch-to-epoch changes (Von Montigny et al. 1995, Ghis- 
ellini et al. 1998). These favor the SSC model in Mrk 421 be- 
cause it is a BL Lac object for which the ratio of thermal (ac- 
cretion disk and broad line region) and synchrotron photons is 
~ . 1 , indicating that the EC mechanism is not important. De- 
tailed modeling of "blue" BL Lacs finds that one-component 
SSC model can generally account for the time-averaged spec- 
tral energy distributions (Ghisellini et al. 1998). Some data 
sets seem to require modifications of the simple model, intro- 
ducing either multiple SSC components or additional external 
seed photons (e.g. Blazejowski et al. 2005). 

We can, however, further decrease the degeneracy among 
proposed physical models by taking advantage of blazars' 
rapid, large-scale time variability with simultaneous X- 
ray/TeV monitoring (e.g. Tavecchio et al. 1998; Maraschi 
et al. 1999; Krawczynski et al. 2001). Different models pro- 
duce emission at a given frequency with particles of different 
energies, cooling times, and cross sections for different pro- 
cesses (e.g. Blumenthal & Gould 1970; Coppi & Blandford 
1990), and thus are in principle distinguishable (Krawczynski, 
Coppi, & Aharonian 2002). For example, the SSC model pre- 
dicts nearly simultaneous variations in both the synchrotron 
and Compton components (see however §4), while other mod- 
els predict more complicated timing (e.g. Ulrich et al. 1997). 

With the possible exception of the X-ray and 7-ray data 
taken on Mrk 501 , the multiwavelength observations on which 
we base our inferences, have often undersampled the intrinsic 
variability timescales (Buckley et al. 1996; Petry et al. 2000; 
Tanihata et al. 2001; Maraschi et al. 1999) and lack a suf- 
ficiently long baseline to make a quantitative assertion about 
the statistical significance of a correlation. Recently, there has 
even been evidence of an "orphan" TeV flare for the Blazar 
IBS 1959H-650 (Krawczynski et al. 2004); a transient 7-ray 
event that was not accompanied by an obvious X-ray flare in 
simultaneous data. 

For what concerns Mrk 421, there have been regular mul- 
tiwavelength campaigns in the last several years, planned 
with an observing strategy focusing on month-long timescales 
(Blazejowski et al. 2005; Rebillot et al. 2006), and in turn 
a relatively sparse time sampling (typically one RossiXTE 
snapshot per night, plus biimed RossiXTE/ ASM light curves). 
These campaigns showed that X-ray and 7-ray brightnesses 
vary "in step" and there is certainly a "loose" correlation, 
and also raised some questions about the need of consider- 
ing additional components to account for the spectra (Blaze- 
jowski et al. 2005), or very high Doppler factors (RebiUot 



et al. 2006). 

The March 2001 campaign remains the experiment with the 
highest density coverage at both X-ray and TeV energies, and 
thus the best dataset to address questions concerning the char- 
acteristics of the variability of the two spectral energy dis- 
tribution (SED) components. Moreover, the brightness state 
achieved during the week of the observations was unprece- 
dented and it remains unparalleled. Preliminary results were 
presented in Jordan et al. (2001); Fossati et al. (2004). 

In this paper we present the summary of the multiwave- 
length observations, with a particular focus on the correlated 
X-ray/TeV variability, also including the TeV data taken by 
the HEGRA (High Energy Gamma Ray Astronomy) tele- 
scope. An account of the HEGRA March 2001 observations 
was published by the HEGRA collaboration (Aharonian et al. 
2002). Giebels et al. (2007) report on additional TeV observa- 
tions with the CAT observatory (Cerenkok Array at Themis), 
simultaneous with the RossiXTE data, not included in this pa- 
per. 

The RossiXTE observations, timing and spectral properties 
are fully presented by Fossati et al. (2008, in preparation; 
hereafter F08). 

The paper is organized as follows: the relevant information 
about the X-ray and TeV observations, and data reduction is 
given in §2. The observational findings are presented in §3. 
and discussed in §4 in the context of the synchrotron self- 
Compton model. Section 4 summarizes the conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. RossiXTE 

RossiXTE observed Mrk 421 in the spring of 2001, dur- 
ing Cycle 6 for an approved exposure time of 350 ks (ObsID 
60145). Observations started on March 18*, 2001, and lasted 
until April V\ 2001, yielding a total of 48 pointings. The 
RossiXTE sampling was very dense until March 25. Dur- 
ing the second week, RossiXTE observed Mrk 421 only dur- 
ing the visibility times for Whipple, whose visibility windows 
were also drasticaUy shortening. In this paper we only present 
and discuss the data obtained between March 18 and March 
25. The journal of this subset of RossiXTE observations is re- 
ported in Table 1 . A complete account of the X-ray campaign 
is presented in F08. 

There are two pointed instruments on-board RossiXTE, the 
Proportional Counter Array (PCA) (Jahoda et al. 1996), and 
the High Energy X-ray Timing Experiment (HEXTE) (Roth- 
schild et al. 1998). 

The PCA consists of a set of 5 co-aligned xenon/methane 
(with an upper propane layer) Proportional Counter Units 
(PCUs) with a total effective area of Ri6500 cm^. The instru- 
ment is sensitive in the energy range from 2 keV to w 100 keV 
(Jahoda et al. 1996). The background spectrum in the PCA 
is modeled by matching the background conditions of the ob- 
servation with those in various model files on the basis of the 
changing orbital and instrumental parameters. 

Since the source was very bright throughout the whole cam- 
paign, we used "bright source" selection criteria^^ to select 
good time intervals (GTl), and all layers of the proportional 
counter units. We did not include the PCUO data^° in this 

OFFSET<0.02, ELY >10, TIME_SINCE_SAA >5. 
In May 2000 the propane layer of PCUO was lost, resulting in a signifi- 
cantly higher background rate and calibration uncertainty. 
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analysis. 

The total exposure time for the PCA was ^263 ks, divided 
over 98 pointing "segments", i.e. intervals with the same 
PCUs on, with individual exposure times ranging between 
144 s and 3.3 ks (we rejected 4 GTIs lasting only 16 s). The 
number of PCUs operational during each pointing varied be- 
tween 1 and 2 (excluding PCUO), with the vast majority of 
cases (97/98) having 2 (for details please refer to F08). 

HEXTE consists of two clusters of four Nal(Tl)/CsI(Na) 
phoswich scintillation counters that are sensitive in the range 
15-250 keV. Its total effective area is «1600 cm^. The 
HEXTE modules are alternatingly pointed every 32 s at 
source and background (offset) positions, to provide a direct 
measurement of the background during the observation, there- 
fore allowing background subtraction with high sensitivity to 
time variations in the particle flux at different positions in the 
spacecraft orbit. Thus, no calculated background model is re- 
quired. The GTIs for the HEXTE data have been prepared 
independently, i.e. not requiring any PCU to be active, thus 
resulting in slightly different on-source times (see Table 1). 
The total HEXTE on-source time was ~282 ks. 

The standard data from both instruments were used for the 
accumulation of spectra and light curves with 16 s time reso- 
lution. Spectra and light curves were extracted with FTOOLS 
v5 . 1 . PCA background models were generated with the tool 
pcabackest, from RossiXTE GOF calibration files for bright 
sources. 

The analysis described in the following is mostly based 
on the 2-15 keV data from the Proportional Counter Array. 
However, the exceptional brightness of Mrk 421 during this 
campaign makes it possible to accumulate a light curve with 
good time resolution also for the higher energy data from the 
HEXTE detector. HEXTE data between 20 keV and 60keV 
were used. 

Light curves for intervals with different active PCUs have 
been rescaled to the same count/s/PCU units by using the rel- 
ative weight for the different effective areas as discussed in 
the RossiXTE GOF website. 

The very high count rate in the RossiXTEfPCA enables 
us to obtain high S/N light curves for several different en- 
ergy bands, and this to study the energy dependence of the 
flux variability and X-ray/TeV relationship. Based on the 
distribution of counts (averaged over the range of observed 
spectral variability), and the relative background contribu- 
tion at different energies, we selected the following four 
bands (expressed in terms of the "hardware" channels) that 
have approximately the same statistics: ch. 5-8 (labeled 2- 
4keV), ch. 10-12 (4-6 keV), ch. 14-18 (6-8 keV), ch. 20-37 
(9-15 keV). Each band comprises on average ^17-24% of the 
PCA counts. For more details please refer to F08. 

We include here some limited discussion involving HEXTE 
data (namely their brightness correlation with the TeV 
dataset). 

A non-variable neutral Hydrogen column density Nh of 
1.61 X 10^" cm~^ (Lockman & Savage 1995) have been used. 
However, since the PCA bandpass starts at about 3 keV, the 
adopted value for Nh does not affect significantly our results. 

2.2. Whipple Observatory 'j-ray data 

The TeV data on Mrk 421 were taken with the Whipple 
10m atmospheric Cherenkov telescope. The Whipple tele- 



scope detects 7-rays by imaging the flashes of atmospheric 
Cherenkov light emitted by gamma-ray induced electromag- 
netic showers. Individual shower images are recorded on a 
photomultiplier tube (PMT) camera, and later analyzed to re- 
ject the cosmic ray background, select 7-ray events, and de- 
termine the point of origin and energy of each detected 7-ray. 
For the period of the 2001 Mrk 421 observations, the cam- 
era consisted of a hexagonal array of 490 PMTs on a 0.12° 
pitch PMT camera (only the inner 379 pixels were used for the 
present analysis) (Finley et al. 2001). Images were character- 
ized by calculating the moments of the angular distribution of 
light registered by the PMT camera. These moments include 
the RMS width and length of the shower images, and orienta- 
tion of the roughly elliptical images. Background cosmic-ray 
showers are rejected using data cuts based on these parame- 
ters according to the procedure described by de la Calle Perez 
et al. (2003). The point of origin of each detected 7-ray can be 
determined from the orientation and elongation of the image 
to with a precision of about 0.12 degree (Lessard et al. 2001). 

The total integrated signal for a shower image (the shower 
size) is roughly proportional to the energy of the primary 7- 
ray. To obtain a better energy estimator, we also correct for 
the impact parameter of the shower. An approximate measure 
of the impact parameter of each shower is obtained by mea- 
suring the parallax angle between the image centroid (point 
of maximum shower development) and the source direction. 
Correcting for the dependence on this centroid distance, one 
obtains an energy estimator for each candidate 7-ray event. 
Monte Carlo simulations of 7-ray showers are used to deter- 
mine the relationship between this energy estimator and the 
actual energy. To determine the energy spectrum for a given 
data set, a histogram of the energy estimator for candidate 7- 
ray events is formed for both the ON-source and OFF-source 
data. The difference in the ON-OFF histograms is then con- 
verted to an energy spectrum by using the effective area func- 
tion calculated by detailed Monte Carlo simulations of the air 
shower and detector. 

For the present analysis, the spectra were reconstructed ei- 
ther using the method of (Mohanty et al. 1998; Krennrich 
et al. 2001) (method A) or the forward folding method of Re- 
biUot et al. (2006) (methodB). In methodA, the combined 
effects of the spectral deconvolution (assuming a log nor- 
mal resolution function) and energy-dependent efficiency are 
taken into account by dividing the binned fluxes by the effec- 
tive area function. The spectral slope is then determined by 
fitting a spectral model to the unfolded spectral data points 
by minimization. Following Krennrich et al. (2001), we 
assume that the Mrk 421 spectrum can be characterized by a 
single power-law with an exponential energy cutoff at a fixed 
value of 4.3 TeV. 

In methodB, instead of fitting to the unfolded spectrum, 
we directly compared the energy estimator distribution of a 
simulated 7-ray data set to the measured distribution for the 
real data. A grid search was used to find the parameters (flux 
normalization and spectral index) that resulted in the best fit 
to the data. To make this computationally tractable, we did 
not repeat the simulations for each trial, but drew our random 
datasets from a weighted simulation database. The measured 
and simulated energy-estimator distributions were compared 
to form a hkehhood for each trial value of the spectral index. 
An adaptive grid search was used to find the best fit value of 
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the flux normalization and spectral index, and to determine 
the confidence interval on the best fit parameter 

For our spectral analysis, it is important to know the energy 
scale to get a reliable value for the flux normalization. For 
this season the trigger condition (which directly impacts the 
energy threshold) consisted of the requirement that 3 adjacent 
PMT signals exceeded a threshold of 32 mV. By comparing 
muon-ring images obtained in the data with those generated 
in Monte Carlo simulations, we obtained a gain caUbration for 
the 2001 season. 

Even in its relatively active state, the 7-ray observations of 
Mrk 421 are often background limited. To adequately char- 
acterize the background of misidentified cosmic -ray events 
we either used an ON-OFF analysis or tracking analysis. For 
these observations, we employed an observing strategy that 
was a compromise between the desire to obtain continuous 
coverage for the time-series analysis, and the need for inter- 
spersed OFF-source control runs to adequately characterize 
systematic variations in background. 

In a typical night, a single ON/OFF pair was taken in which 
a 28-minute ON-run run with the 10m telescope tracking 
Mrk 421 was followed by another 28-minute observation that 
tracks the same range of azimuth and elevation angles, but off- 
set by 30 minutes in right ascension. This mode of data taking 
reduces systematic errors from difference in sky brightness, 
atmospheric conditions and instrument variations but reduces 
the duty cycle of observations. Other runs taken in the same 
night were typically acquired in a staring or tracking mode 
where the telescope continuously tracks the source without 
interruption for off-source observation. 

Data selection cuts based on the RMS width and length of 
the shower images were used to determine candidate 7-ray 
events. The angle between the major axis of the Cherenkov 
image and the line connecting the image centroid to the angu- 
lar position of the source (designated a), is used to define the 
signal and background region in the image parameter space. 
In each ON-source run, the events whose major axis points 
to the position of the source in the field of view are used 
to estimate the signal, and those misaligned events pointing 
away from the center are used to determine the background 
de la Calle Perez et al. (2003). The number of candidate 7- 
ray events was determined by subtracting the number of on- 
source events from off-source events. Significances were de- 
termined by the method of Li & Ma (1983). 

Since the data were taken over a range of zenith angles, a 
first order correction was made to take into account variations 
in energy threshold and effective collection area with zenith 
angle. For the light curves presented in this paper we choose 
the simple empirical method of normalizing to the flux from 
the Crab Nebula at the corresponding zenith angle. Like rel- 
ative photometry, this method cancels systematic errors and 
should provide a better method for eventually combining data 
with data obtained by different detectors calibrated by differ- 
ent Monte Carlo simulations. However, this method results 
in a systematic (second order) error if the Mrk 421 spectrum 
differs significantly from that of the Crab Nebula. 

The Crab Nebula can be used to empirically quantify the 
combined effect of these changes for a source with the same 
spectrum A^(£') - E'-^ ^'' (Hillas et al. 1998). Fi sure 1 shows a 
fit to the Crab rate as a function of zenith angle for 50 ON/OFF 
runs from October 25, 2000 to the end of the 2001 observing 
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Fig. 1 . — Dependence of the observed rate of gamma rays from the 
Crab Nebula as a function of zenith angle, as observed in the 2000- 
2001 season. The fit function indicated by the dashed curve was used 
to normalize the raw Mrk 421 rates. 



season. 

The average spectral index of Mrk 421 has been measured 
to have a similar value to the Crab nebula in the Whipple 
energy range, but with evidence for some spectral variabil- 
ity (Piron et al. 2001; Krennrich et al. 2002; Aharonian et al. 
2002). Since the energy threshold only varies by a factor of 
roughly 30% over the range of zenith angles of our observa- 
tions, we estimate a systematic error in the flux of roughly 
10% if we assume that the spectral index varies by AT ^ 0.3. 
But we lack the ability to determine the spectral index to this 
precision for most individual runs, and thus subsequently ad- 
mit the possibiUty of some systematic errors, and normalize 
our observations to a functional form for the zenith-angle de- 
pendent rate from the Crab Nebula. 



7.423 secp-iO) 



[exp(0.5sec(6'))sec2(6i))] 



1.49 



(1) 



This function is derived from an empirical fit to data taken 
on the Crab Nebula at various zenith angles during the same 
observing season (see Figure 1). 

Total observations for the March 2001 campaign came to 
10.3 hours of ON/OFF data and 37.7 hours of TRACKING 
data for a total of 48 hours of data, only a fraction of which 
are covered in this paper. 

For our multiwavelength studies we chose three arbitrary, a- 
priori bin widths i) one corresponding to a single day of obser- 
vations, ii) one corresponding to a single 28 minutes data run, 
and iii) a shorter bin width of ?a4 minutes. The latter choice is 
arbitrary, but was derived from past analyses of strong flares 
as giving a good compromise between statistics and temporal 
resolution. 

2.2.1. HEGRA Observatory ^-ray data 

The HEGRA TeV data, light curve and spectra, utilized in 
this paper have been previously published by Aharonian et al. 
(2002). Please refer to the original paper for details concern- 
ing the data reduction. 
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2.3. Combining Whipple and HEGRA j-ray data 

The locations of the Whipple and HEGRA telescopes, 
in Arizona and the Canary islands, respectively, separated 
by approximately 6 hours, make it possible to achieve un- 
interrupted coverage of Mrk 421 during the spring, when the 
target is observable at "small" zenith angles for up to 7 hours 
each night from each site. The visibility windows thus com- 
plement each other very well (see Figures 2 and 3), but don't 
overlap. For this campaign we achieved an unprecedented 
coverage of the target: as illustrated in Tables 2 and 3 the total 
net exposure time for the two TeV observatories was of about 
62 hours, over seven nights. Ignoring gaps shorter than 1.5 
hours, the on-source time was of about 72 hours, i.e. more 
than 10 hours per night. 

The combination of data from different instruments is al- 
ways a very delicate step, most robustly addressed by compar- 
ing data taken simultaneously, which is not possible for Whip- 
ple and HEGRA. The next best option is to cross-calibrate the 
data using a standard candle as reference, typically the Crab 
Nebula for high energy emission. The main issue concerning 
the combination of the data from the Whipple and HEGRA 
telescopes arises from the fact that i) they gather data with 
different lower energy threshold, 0.4 and 1 .0 TeV respectively, 
and that ii) Mrk 421 TeV spectrum is in general (significantly) 
harder than the Crab's. To illustrate the problem, let's take 
spectral indices F = 2.5 for the Crab, and F = 2.2 for Mrk 421 
(Krennrich et al. 2002; Aharonian et al. 2002). Ignoring the 
fact that a detector response is energy dependent we can com- 
pare the ratio between the Mrk 421 and Crab fluxes computed 
for different energy thresholds, namely 0.4 and 1 TeV. The 
flux above an energy E for a simple power law with spectral 
index F is: 

F(£,F) = ^Ml^; (2) 



Therefore 



FmrkiE) ■■ 



F(E,T„,k) No. 



F-1 



,mrk ^ Crab ' 



1 



F(E , Tcrab) Nofirab ^mrk " 1 



(3) 



and comparing data taken at two different thresholds yields. 



/ = 



FmrkiEB) 



(4) 



For £^ = 1 .0 TeV (HEGRA) and £b = 0.4 TeV (Whipple), this 
ratio is 1.32, i.e. a spectrum yielding a Whipple flux of 1 Crab 
wiU be observed at 1 .32 Crab by HEGRA. Moreover, this ef- 
fect is a function of the spectral indices. While we can safely 
consider the Crab spectrum non variable, significant variabil- 
ity is observed in Mrk 421 (Krennrich et al. 2002; Aharonian 
et al. 2002). The value of / = 1 .32 for die HEGRAAVhipple 
flux ratio obtained for T„rk = 2.2 becomes / =1 .44, 1 . 10, 0.83 
for Tmrk = 2.1,2.4,2.7, respectively. Hence, HEGRA data 
would also show a larger variance. Finally, it is worth noting 
fliat the ratio between the March 2001 Whipple and HEGRA 
mean fluxes, and between their variances, are consistent with 
what expected based on the simple arguments just illustrated 
and the observed spectral indices. 

Since it is not possible to derive from "first principles" a 
robust way to convert Whipple and HEGRA light curve mea- 
surements into each other, we conclude that the most robust 
procedure at hand is to scale the data of one telescope to the 



same mean and variance of those of the other one. We adopt 
Whipple as reference and scale the HEGRA fluxes. We deem 
that this approach is rehable, primarily because of the large 
size of the two datasets: they both span seven days (with al- 
ternating visibiUty windows), and have comparable variability 
amplitudes, thus both should be sampling an equally repre- 
sentative subset of the source phenomenology. Being a linear 
conversion, this procedure can not correct for the slight non- 
hnearity of the Whipple vs. HEGRA flux relationship intro- 
duced by the brightness-hardness correlation of the spectra. 
In order to mitigate the effect of outliers, and because of the 
intrinsic exponential nature of the source variability, we per- 
formed the seating of mean and variance on the logarithm of 
the Ught curves. 

2.4. Optical data 

UBVRI optical monitoring was performed using the 
Harvard-Smithsonian 48" telescope on Mt. Hopkins. Data 
were analyzed using relative aperture photometry, using com- 
parison star #3 as listed in Villata et al. (1998) Galaxy 
background light was subttacted using the simple empiri- 
cal method described in Nilsson et al. (1999). Here we use 
Nilsson et al. (1999)'s determination of the contribution of 
Mrk 421 galaxy background light in the R band, but extrapo- 
late this to the V band using the - V color of the host galaxy 
as given in Hickson et al. (1982). The errors shown on the 
Ught curve are the systematic uncertainties, which dominate 
the small statistical errors. These errorbars were determined 
from the measured variance of the reference stars with re- 
spect to each other, and may not include other effects such 
as the bleeding of starhght from the bright stars SAO 62387 
and SAO 62392 tiiat lie about 2 arcmin from Mrk 421, in tiie 
48" telescope field of view. Other systematic effects come 
from the relatively high level of galaxy light from Mrk 421 
and the very nearby satellite elliptical galaxy. These prob- 
lems, coupled with the lack of good reference stars combine 
to make Mrk 421 one of the more difficult BL Lac objects for 
optical monitoring; we estimate that this measurement should 
be given a systematic flux uncertainty of about 15%. 

3. TEMPORAL ANALYSIS 

3.1. RossiXTE/Whipple+HEGRA data overlap 

The principal statistics regarding the quaUty of the X-ray/7- 
ray overlap are summarized in Table 2 and 3. There we re- 
port two different measures of the overlap between the two 
telescopes. The fractions listed in columns 5 and 8 refer to 
the absolute "covering factor" by RossiXTE of the on-source 
time of Whipple or HEGRA. These numbers give an idea of 
how "representative" the collected X-ray data are of the ac- 
tual X-ray brightness. In columns 6 (and 9), instead, we re- 
port the fraction of Whipple (and HEGRA) data points (from 
the run-by-run light curves) that have a non-zero overlap with 
any fraction of the RossiXTE/PCA GTIs: the actual exposed 
fractions for the individual runs range between 7% and 100% 
(with average of the order of the value reported in column 5 
and 8). Although the RossiXTE schedule was optimized to 
observe Mrk 421 as much as possible^', we still missed about 

2' Between Match 19 and 25 RossiXTE observed Mrk 421 for 94 out of 
102 orbits, and in only one case there were two consecutive orbits missed. 
One missed orbit corresponds to a ~ 8.5 ks gap. 
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Fig. 2. — Simultaneous optical (V band, bottom), X-ray (2-10 keV, middle) and TeV 7-ray {E > 0.4 TeV, but see §2.3, top) (top) light curves 
for Mrk 421 for the March 18-March 25 period. RossiXTEfPCA data are shown here in 256 s bins. HEGRA data (dark triangles) are integrated 
over 1800 s bins, Whipple data (white circles) over 1680 s bins. HEGRA data precede Whipple's. The optical data have been rebinned to yield 
a s/n ratio of at least 8, but with bin length not exceeding 1500 s. The logarithmic scales span a factor of vTo, 10 and 100 for the optical. 
X-ray and 7-ray light curves, respectively. The length of axes scale accordingly ( x 2 between them), so that relative amplitude variability can 
be directly compared. 



1/4 of the Whipple runs, and during each Whipple/HEGRA 
« 0.5 hr integration window X-ray data were collected only 
for about 50% of the time. Since the X-ray and 7-ray bright- 
ness can vary significantly on timescales faster than w 0.5 hr, 
we may expect this to affect the X-ray/TeV flux correlation. 
We carefully inspected the X-ray light curves to assess the im- 
pact of the un-even coverage and concluded that the effect is 
not important. The largest X-ray variation during any 0.5 hr 
interval is ^30%, with an average of ^15%. 

3.2. Light curves 

The resulting light curves for the full week are shown in 
Figure 2. The overall brightness ranges are approximately 
xO.5, xlO, xlO for optical. X-ray and 7-ray respectively. 
Although the time coverage was the best possible and the 
datasets encompass several days and multiple flares, these 
ranges can be affected by the lack of continuous coverage by 
the ground based telescopes (e.g. there are no optical or 7-ray 
data for the time of the lowest X-ray flux, on March 19). 

Although the large gaps in the 7-ray light curves may limit 



the interpretation, some effects are obvious. The source shows 
strong variations at both bands, and these variations are highly 
coiTclated, with features in one band generally showing up in 
the other band as well. This is clearly illustrated in Figure 3, 
where we plot 50 ks sections of each night, centered on the 
time window covered by the TeV observations. 

The best example is the March 19 observation, when a well 
defined, isolated, flare was observed both in the X-ray and 
7-ray bands, from its onset through its peak (more in Sec- 
tions 3.3 and 3.4). In several other nights, despite the lack of 
other major "isolated" features in the light curve(s) that can 
be provide an unambiguous reference for the comparison, the 
variations at X-ray and 7-ray are clearly correlated. Given the 
high degree of variability (rapid and high amplitude in both 
bands) it is always possible that intrinsically un-correlated 
flares in the two bands end up being simultaneous by chance. 
However, the unprecedented level of detail (i.e. time resolu- 
tion) and extension of the simultaneous coverage, allows us to 
"match" several relatively minor features in the light curves, 
and effectively confirm the coiTclation of the variations in the 
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Fig. 3. — Simultaneous 2-10 keV X-ray and TeV (see text) 7-ray light curves for individual nights. Gray triangles are HEGRA data, in « 1800 s 
bins. White circles are Whipple data, integrated over ~1680 s bins. Dense dark points are RossiXTEfPCA, in 128 s bins. The shaded boxes 
represent the average and variance of the X-ray data for each (longer duration) TeV bin, that are the values used in the Flux-Flux correlation 
analyses. The rate scales for the X-ray data are on the left Y-axes, and the flux scales for the TeV data on the right Y-axes. The time span is the 
same for all panels, 50 ks. The vertical scales are not the same in all panels, but are adjusted to show each day in the best possible detail. The 
X-ray dynamic ranges are (time ordered) x6, x5, x2, x4, x4, x4, x4. In order to allow for an easier comparison of the relative variability 
amplitude, in all panels, the Y-axis range for the 7-ray light curve is the square of that used to plot the X-ray data. The source shows strong, 
highly-correlated variability in both energy bands, with no evidence for any interband lag (see however §3.3). 



X-ray and 7-ray bands. 

In particular, within the detail allowed by the coarser TeV 
sampling, the X-ray and 7-ray light curves seem to track each 
other closely in all cases when the sampling is good, namely 
for the nights #1, 2, 4, 5, 7. For nights #3 and 6 there seem 
to be significant deviations from this general trend. It is, how- 
ever worth noting that these two nights correspond to the high- 



est brightness level in the X-ray, and the X-ray light curves 
show several fast variations (intra-orbit), that may not have 
been sampled properly by the TeV observations (e.g. in #3), 
or may in fact not have a TeV counterpart (e.g. #6, where the 
TeV data seem to have good signal to noise). In a broad sense, 
even in #3 and #6, the data are consistent with correlated vari- 
ations in the two energy bands. 
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Fig. 3. — Continued. 



In the next two subsections (§3.3, §3.4), we are going to 
examine in detail the properties of the X-ray/TeV correlation 
from two complementary point of views, phase and ampli- 
tude. 

3.3. X-rayfTeV interband lags 

Cross-correlation functions were measured to quantify the 
degree of correlation and phase differences (lags) between 
variations in the X-ray and 7-ray bands, using the discrete 
correlation function (DCF) of Edelson & Krolik (1988). For 
the X-rays we consider two different bands (2-4 keV and 
9-15 keV) at the usable ends of the PCA bandpass, to ex- 
plore the possible energy dependence of the correlation. The 
results for the whole-week and March 19 flare data are shown 



in Figures 4a-d. It is worth noting that, as discussed in F08, 
variations in the harder X-ray emission seem overall to lag 
those in the softer PCA band. 

3.3.1. Full week-long dataset 

For what concerns the full week dataset (Fig. 4a,b), we do 
not find a measurable lag, with either X-ray band. The statis- 
tical properties of the DCF are complicated by the presence 
of several regular patterns in the data trains (chiefly the diur- 
nal gaps in the ground based TeV data, and the orbital gaps 
of the /?o55/XTE data). The DCFs peak at zero-lag, and the 
correlation coefficients are quite high, ~ 0.8, despite the poor 
sampling (the combination of the large diurnal gaps due to 
the ground based visibility, and the 7?oii/XTE orbital gaps 
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yields an efficiency of about 1/2 x 3/4 2± 40%). By means 
of simple data-based simulations comprising flux randomiza- 
tion (e.g. Peterson et al. 1998), and the effect of introducing 
a shift in one of the timeseries, we estimate an upper limit of 
|r| < 3 ks on the value of the soft-X-ray/7-ray lag, possibly 
smaller for the harder-X-ray/7-ray case. 

3.3.2. The March 19 flare 

Beside the full dataset, we focused our attention on the iso- 
lated outburst of March 18/19 that uniquely comprises many 
favorable observational characteristics, namely i) the best TeV 
coverage, ii) the least RossiXIE data gaps, and iii) the best 
/?oiijXTEAVhipple overlap, iv) the largest brightness excur- 
sion (in both bands, x 10 in 7-rays, x3 in X-rays). Because 
of this indeed rare combination of properties, the March 18/19 
flare's DCF is probably not significantly affected by the sam- 
pUng. Because of its relative isolation from other outburst, 
and its large amplitude, we can regard this event as a rather 
clean "experiment", providing us a good view of the variabil- 
ity mechanism at work. On the other hand, results obtained 
for this particular event may not necessarily be representative 
of all flares. 

In this case we used the Whipple data in their shortest avail- 
able time binning, 256 s bins (an example of this is the peak 
region of the March 19 burst shown in Figure 5c). The DCFs 
for the two different X-ray bands, computed over 1024 s time 
steps, are plotted in Figures 4c and 4d. The DCFs for both 
X-ray bands show a high correlation coefficient, peaking at 
0.84 (at a lag T « +2 ks) for the soft X-rays, and at 0.88 (at 
zero lag) for the harder X-rays. 

The most remarkable feature is that there seems to be a hint 
for the 7-rays lagging the softer X-rays, while being "syn- 
chronized" with the harder X-ray photons. A thorough anal- 
ysis and characterization of the properties of the X-ray vari- 
ability is discussed in F08. In Figure 5a,c, we just show a 
20 ks section of the three light curves for March 19, centered 
around the possible peak position. Although the peak of the 
outburst was not directly observed in X-rays, we can make the 
following heuristic arguments concerning the possibility, and 
value, of a interband lag between the softer X-rays and the 
TeV data. 

3.3.3. Constraining the lag for the March 19 flare 

By exploiting the statistical knowledge of the characteris- 
tics of the X-ray variability (see F08), during this campaign, 
we can try to assess the probabiUty that a lag at the flare peak 
in fact exists between the TeV and the soft X-ray light curves. 
The idea is to assign a probability distribution for the soft X- 
ray peak to occur at times fpeak during the data gap, and then 
use it to evaluate the probability that the peak in fact occurred 
within the time interval comprising the peak of the TeV out- 
burst, fliat is at r - Tref > 30.25 ks. 

The basic building blocks for these probability estimates 
{^(?peak)} are the observed distributions of doubling and halv- 
ing times, V{t2) and P(ti/2), derived from the entire week 
long soft X-ray dataset. With these we assign the probabil- 
ity for the peak of the X-ray Ught curve to occur at a certain 
time and brightness (fpeak,fpeak) within the data gap, by tak- 
ing the joint probability of having the T2 and T1/2 required to 
reach each trial (fpeak,^peak) position "moving" from the left 
(i.e. before, ?bpj^bp) and right boundaries of the data gap (i.e. 



after, fap,Fap). 

'^(^peak, ^iieak) ~ 'P(T2(?bp , ^bp ; ^peak, F^skf) ■'Pin/l {tap ) ^ap '1 ^peak ) ^pe 

(5) 

Since we are not interested on Fpeak this distribution is then 
summed over all Fpeak to yield just ^(fpeak)- 

The P(t2) and P{t\/2) adopted in this analysis were de- 
rived from the doubling and halving times from all data pairs 
whose separation in time ATij is between 0.25 and 5 ks. We 
restricted our samphng to this subset of data pairs because we 
wanted the distribution to be representative of the same type 
of variations that could have occurred during the data gap, 
which spans ~3 ks. The inclusion of larger pair separations 
would "spuriously" bias the probabiUty towards large values 
of T2 and Ti/2- On the other hand, relaxing the limit on the 
minimum time separation picks up very fast variations, which 
are not relevant for this analysis, because their influence on 
where the peak could fall is marginal (given their limited am- 
pUtude), and their overall contribution is already taken into 
account (smoothed out) by the "slopes" measured on longer 
timescales. 

We have performed this analysis with different choices of i) 
the allowed range of A7ij, and of ii) the "starting" points on 
both sides of the gap (namely we checked points up to ±2 ks 
from the gap). The results do not change significantly. 

The probability distribution for fpeak resulting from this 
analysis is shown in Figure 5c. The average of the all the dif- 
ferent tests yields a probability of ^(Jpeak > 30.25) ~ 1 -2% 
for the flare peak to occur later then T-T^^f = 30.25 ks, i.e. 
within the Whipple peak time interval. The most probable 
fpeak estimated by this method is fpeak = 28.6 ± 0.8 ks (Icr), 
two sigma below the "first possible time" for the peak of the 
TeV flare. 

We can push this type of analysis a little further to esti- 
mate the most likely value for the lag between soft X-rays 
and TeV. In order to do this we need to assign a probability for 
the time of the TeV peak 7'(fpeak.Tev)- We tried the following 
simple distributions for 7'(fpeak,Tev): i) a uniform distribution 
within the 1680 s integration window, ii) a "tent" function 
centered on the top interval and going to zero at its bound- 
aries, iii) a "tent" function centered on the top interval, but 
extending over half of each of the two adjacent intervals (i.e. 
, T - Tref ~ 29.7 - 3 1 .7 ks). The convolution of the 7^(fpeak,x) 
with 7'(fpeak,Tev) shifted by r gives the probabihty for a given 
lag T. The result does not change significantly with the differ- 
ent choices i)-iii), and it is r = 2.06;!;q ^ ks (Icr). 

It is important to stress that this analysis rests on a few as- 
sumptions, that we deem reasonable, that are here summa- 
rized. 

• The statistical properties of the X-ray variability change 
on a timescale longer than our experiment. In this respect 
we checked that the distribution of T2 and T1/2 for different 
subsets of the week-long dataset are consistent with each 
other. 

• The power spectrum of the variations in X-rays and TeV is 
such that the there is only a negligible probability that the 
peak of the X-ray fight curve occurred before or after of 
the data gap, and that of the 7-ray light curve during one 
of the earlier or later integration windows. For instance we 
rule out that the TeV flare peak could have been reached by 
means of a very fast and very large amplitude variation (a 
spike not resolved, and smoothed out, by the coarse Whip- 
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Fig. 4. — Cross correlation between the X-ray and the TeV hght curves, (a) 2-4 keV vs. TeV (Whipple+HEGRA) for the whole campaign 
(computed over 2048 s bins, from X-ray data on 256 s bins, and TeV data on ~750-900 s bins), (b) 9-15 keV vs. TeV (Whipple-FHEGRA) 
for the whole campaign (computed over 2048 s bins, from X-ray data on 256 s bins, and TeV data on ~750-900 s bins), (c) 2^ keV vs. TeV 
(Whipple) for the night of March 18-19 (the flare of Figure 3a) (computed over 1024 s bins, from X-ray data on 128 s bins, and Whipple data 
on 256 s bins), (d) 9-15 keV vs. TeV (Whipple) for the night of March 18-19 (computed over 1024 s bins, from X-ray data on 128 s bins, and 
Whipple data on 256 s bins). 
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Fig. 5. — a) March 19 light curves for Whipple 7-ray (connected white circles, ~ 1000 or 1680 s bins), and two X-ray bands (dark squares 
for 2-4 keV, gray triangles for 9-15 keV, both in 256 s bins). The rate scale for the X-ray data is on the left Y-axis, and the flux scale for the 
Whipple data on the right Y-axis. To allow for an easier comparison of the relative variability amplitude, the Y-axis range for the 7-ray light 
curve (x 16) is the square of that used to plot the X-ray data (x4). b) Same as a) panel, except that here the TeV (Whipple) light curve data 
have been adaptively rebinned from the 256 s data to a signal-to-noise ratio of at least 6. c) Probability distribution of soft X-ray flare peak time 
derived from the general statistical properties of the short term variability (see text, §3.3). The vertical dashed line marks the leftmost boundary 
of the Whipple time interval comprising the TeV flare peak. 
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Fig. 6. — Plot of the Whipple+HEGRA TeV flux vs. the X-ray count rate in different energy bands (see labels), for the entire week. Different 
gray shades correspond to different nights. For reference, in each panel are shown segments indicating different slopes for the relationship 
between the plotted fluxes. 



pie binning), during one of the two Whipple bins falling 
during the gap in the RossiXTE data. Moreover, higher 
sampling Whipple light curves (see Figure 5b) provide a 
further constraint on the probability, and characteristics, 
of this type of extreme event. For what concerns the X- 
rays, we have the possibility of investigating in more detail 
the properties of the variability on fast(er) timescales (see 
F08). A broad assessment of the reliability of our assump- 
tion can be made by considering the likelihood of a large 
amplitude variation on a timescale shorter than e.g. 250 s, 
the cut-off we applied to our sampling of T2 and ti/2- The 
analysis of the fractional rate variability for Af between 
32-250 s shows that the probability for a AF/F > 20% 
is only ^ 6%. 

• We would also like to point out that it would be desirable 
to use not simply the probability distribution for the r's, 
but the probability for a given change in rate AF/F for 
each given t. However, despite the size of the RossiXTE 
dataset, it is not possible to have a good enough sampling 
for V(—,t), to constitute a significant improvement over 
the uncertainty inherent in the assumption that all AF/F 
are equally probable for a given r. 

The same analysis performed for the 9-15 keV light curve 
yields a fpeak = 30.0 ±0.7 ks (Icr), a 7'(fpeak > 30.25 ks) ~ 39%, 



and an estimate of the lag of the TeV peak of r = 0.73 ± 
0.80 ks, i.e. no measurable lag. 

3.4. X-ray vs. TeV flux correlation 

Comparison of the variability amplitudes (as opposed to 
phases) offers different constraints. As clearly illustrated in 
Figures 2 and 3, the source shows stronger variability in the 7- 
rays than in the X-rays: in fact, in all panels the flux scale for 
the TeV data spans a range that is the square of that of the rate 
scale used for the RossiXTE/PCA data, and the light curves 
run in parallel. This is confirmed in Figure 6 that shows 7-ray 
flux as a function of X-ray count rate in different X-ray energy 
bands. The TeV data are binned on approximately 28-minute 
runs. The RossiXTE count rates correspond to the average 
over intervals overlapping with the TeV observations (as in- 
dicated by the shaded boxes in Figure 3), and the error bars 
represent their variance (height of the shaded boxes). 

We fit the log-log data with a linear relationship (i.e. F^ oc 
F^), which provides a satisfactory description in all cases. 
The best fit slopes are reported in Table 4 (top row), along 
with their errors. We also analyzed the X-ray/TeV correla- 
tion for different sections of the campaign, and individually 
for each night, with the intent of looking for possible varia- 
tions. Individual nights plots are shown in Figure 7. In Ta- 
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Fig. 7. — Plot of the Whipple+HEGRA TeV flux vi. X-ray 2-10 keV count rate for each individual observation night, and the combination of 
nights 1+2 and 5 + 6 + 7 . Axes range is x 1 in all panels except for those involving the March 1 9 (night 1 ) data, whose variation range is larger 
(x30). Some of the symbols are shaded in gray to help recognize different sets. 



ble 4 we report the best fit correlation slopes for the best sin- 
gle nights, and for a few combinations of consecutive nights 
(#1-1-2, #4-1-5, #6-1-7, #5-1-6-1-7). 

It is worth noting that again the March 18/19 (day 1) X- 
ray/7-ray observations provide the best case study, for the 
large amplitude of variability, likely ensuring us that the ob- 
served amplitude is close to the intrinsic one. The presence, 
and contribution, of a steady (variable on longer timescale) 
emission diluting the flaring one could alter the perceived am- 
plitude of flares. This is a long-standing issue that is difficult 



to address, but in this respect the March 18/19 flare is a unique 
event. 

The unprecedented quality of this dataset enables us not 
only to establish the existence of the correlation between the 
TeV and X-ray luminosities, but also to start unveiling some 
of its more detailed characteristics, e.g. its evolution with 
time. The emerging picture is complex. There are several 
observational findings that we would like to point out. 

• The first, most direct and general, observation is that the 
TeV flux shows a definitive correlation with the X-ray rate. 
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Fig. 8. — a): X-ray vs. 7-ray one-day averaged brightnesses. The error bars represent the variance during the interval, which can be 
considerable. The correlation is approximately linear (see Table 4). Numbers refer to the campaign night sequence, b): boxes (approximately) 
representing the regions of the diagram occupied by the data of each individual night. The combination of steep(er) intranight and flat(ter) 
longer term, due to shift of the "barycenters", correlations is more easily shown. Both plots are on the same axes scale and range. Box contours 
are solid for nights 1, 2 and 7; long-dashed for 3 and 4; dotted for 5 and 6. 



for all X-ray energy bands (see Figure 6). Considering the 
entire week-long dataset, 105 data pairs, the correlation is 
approximately linear (see Table 4). The same is apparent 
when looking at the data binned over 1-day timescale. Fig- 
ure 8a. 

• A more careful inspection of the flux-flux diagrams sug- 
gests however a richer phenomenology. In fact, we may be 
observing a series of parallel "flux-flux paths", individu- 
ally obeying a steep (e.g. quadratic) trend, but that taken 
together produce a rather flat envelope producing the lin- 
ear trend emerging for the global cases, because of a drift 
of their barycenters. There is indeed a secular increase of 
the source brightness over the course of the campaign, and 
it seems to be more enhanced in X-ray. Its amplitude is 
of the order of intranight brightness variations, thus alter- 
ing the X-ray/7-ray coiTelation on longer timescales. Fig- 
ure 8b shows how the regions covered by nightly data shift 
from day to day, while broadly maintaining an approxi- 
mately quadratic intranight flux correlation trend in most 
cases. 

There is thus an intriguing hint that there might be a 
split between the correlation observed on short (hours) 
timescales and that apparent on longer (days) timescales, 
once faster variations are smoothed out. 

• There may be two different (luminosity related) regimes 
for the X-ray/TeV flux correlation. By splitting the data 
in two sections of significantly different average bright- 
ness level, days 1+2 (with or without the pre-flare noisy 
HEGRA data section), and days 6+7 (or 3+6+7), we note 
that that source seems to exhibit two different behaviors: 
the TeV vs. X-ray relationship is significantly steeper for 
the day- 1+2 subset, with values of 77 for all 4 PCA en- 



ergy bands larger than 77 = 1.82(±0.12), versus all values 
smaller than rj = 1 .03(±0. 14) for days 6+7 (Table 4). 

• For two nights (1 and 5) the flux-flux diagram is very tight, 
with all points lying on a very narrow path. In these cases 
the TeV flux increases more than linearly with respect to 
the X-ray rate. For the flare of March 19 the coiTelation is 
"super-quadratic" at all energies (Table 4). Moreover, for 
these two nights the light curves encompass a full flaring 
cycle, i.e. we can follow the complete evolution of an out- 
burst, rising and decaying. In both cases the paths of the 
rising and decaying phases in the flux-flux diagram over- 
lap perfectly. 

• There is no significant change of the slope of the correla- 
tion with the choice of X-ray energy band, except for the 
case of the full-week dataset. A flatter coiTelation slope 
for harder X-rays would be expected because of the in- 
trinsically higher amplitude of the variability of the syn- 
chrotron component towards higher energies (if we are al- 
ready above the peak energy) (e.g., Fossati et al. 2000a). In 
fact, the relative variance, af/{F), increases with energy, 
changing from ~ 0.45, to 0.48, 0.52, 0.56 for 2-4, 4-6, 
6-8, 9-15 keV respectively. This change fully accounts 
for the flattening of the X-ray/7-ray correlation slope. The 
effect is not observed for smaller subsets of data probably 
because of the lower statistics. 

The departure of the 20-60 keV band from this trend could 
instead be justified by considering that the flux in this band 
may comprise a contribution from the onset of the inverse 
Compton, which could be regarded as constant because it 
would be varying on much longer timescales. However, 
this hypothesis does not seem to be supported by the data, 
because, though very noisy and with limited energy lever- 
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age, the HEXTE data are consistent with the extrapolation 
of the steep PCA power law. Alternatively it is possible 
that the difficult background subtraction of the low-count- 
rate HEXTE data reduces the intrinsic dynamic range of 
the X-ray variations, thus steepening the correlation. 

These observational findings have important implications 
for the physical conditions and processes responsible for the 
variability in the scattering region, as discussed in §4. 

3.4.1. Comments 

Before we proceed to discuss the observational findings, we 
would Uke to put forward a few additional comments concern- 
ing some aspects of the derivation and interpretation of the 
flux-flux correlation. 

• For simplicity we performed the brightnesses correlation 
analysis using count rates for the RossiXTE data. A proper 
conversion to flux units requires to fit a model to the data 
for each short sub-interval, and it would introduce a dif- 
ferent source of uncertainty. We tested the correspondence 
between count rates in RossiXTE/PCA bands and model 
fluxes for a broken power model, with different spectral 
indices, and break energy positions, covering the range of 
values observed in March 2001 (for full account of the 
spectral analysis please refer to F08). For the 2-lOkeV 
band, the correlation between count rate and flux is slightly 
tilted, in the sense of slightly less than linear increase of 
the flux with rate, Flux ^ Rate" '*. This would thus further 
steepen the TeV/X-ray flux-flux correlation if computed 
with the X-ray flux. The effect is small and it is not present 
when narrower energy bands are considered. 

• Rebinning the data alters the variance of the light curves, 
and if the effect is different for X-ray and 7-ray (namely if 
their intrinsic power spectra are different), it could poten- 
tially bias the slope of the correlation. The comparison of 
the change of variance of X-ray and 7-ray (starting from 
the 256 s-binned one when possible) light curves for dif- 
ferent rebinnings, suggests that the effect is at most of the 
order of 10%. The effect is small in comparison with the 
overall range spanned by the data, which is of the order of 
a factor of at least five for the week long dataset. Hence, 
we deem the effect of the choice of time binning on the 
determination of the flux-flux correlation slope not signifi- 
cant. 

• Since we are measuring the fluxes in limited energy bands, 
the slope of the relation depends also on the position of the 
synchrotron and 7-ray peaks with respect to the observed 
energy bands. The reason is that as the peak moves from 
lower frequency into the bandpass of a detector, a small 
change in the peak position yields a larger variation of the 
flux. A simple shift in frequency would be degenerate with 
a true increase in luminosity. Once the spectral peak falls 
within the bandpass, and it is shifting within it, this "spu- 
rious" effect becomes un-important. In a very simplified 
case, taking MrkSOl as test SED, Tavecchio et al. (2001) 
showed that the 7-ray vs. X-ray flux relationship predicted 
for variability simply due to a change in the maximum par- 
ticle energy (and in turn synchrotron peak energy), can 
vary between flatter-than-linear to steeper-than-quadratic 
(the effect was however enhanced by the fact that the au- 
thors compared monochromatic fluxes). Katarzynski et al. 



(2005) performed a thorough analysis of the effect of the 
position of observed energy bands with respect to the syn- 
chrotron or IC peak energies in the context of the X-ray vs. 
7-ray brightness correlation, and found that it can change 
the slope over a broad range of values, including linear and 
quadratic. This apparent freedom is however lost if data 
following the full evolution of a flare are available. In fact 
their conclusion with respect to an outburst developing like 
that of March 19 is that explaining the observed correlation 
by means of specific choices of spectral bands is problem- 
atic and it would require very contrived assumptions. 

The characteristics of the X-ray variability itself seem to 
evolve during the campaign. In particular, it is important to 
recall that the spectrum becomes significantly harder over the 
course of the week-long campaign, accompanying a gradual 
brightness increase. Rather than a caveat this is probably a 
point in support of the apparent change of X-ray/7-ray behav- 
ior between the first and second part of the week. The spectral 
analysis of the March 2001 Whipple data, reported separately 
by Krennrich et al. (2003), showed that the TeV spectra also 
significantly hardened between March 19 and 25. The spec- 
tral indices for a power law fit with exponential cutoff (fixed at 
4.3 TeV) shift from T ~ 2.3 to T ~ 1.8 (±0.15), i.e. suggest- 
ing that the IC peak moved from below to within the Whipple 
bandpass (i.e. in the latter case the 7-ray emission would peak 
at about 1 TeV). RossiXTE spectra present a similar picture of 
the X-ray evolution, namely that the synchrotron peak shifted 
into the PCA bandpass. Broken power law fits show that the 
lower energy spectral index becomes harder than F =2 (FOB). 
Unfortunately even with the available statistics, because of the 
limited energy leverage, it is not possible to pinpoint robustly 
the energy of the synchrotron peak and its evolution (as was 
the case with BeppoSAX). 

If the peak of one component (synchrotron or IC) moves 
into the observed band, we would then be observing the vari- 
ations of a lower, possibly below peak, section of the electron 
spectrum, instead of the more highly variable higher energy 
end. Depending on whether this happens to both peaks or just 
one, we expect to observe a different phenomenology: e.g. if 
this happens only for the TeV band the correlation with the 
X-ray data should become flatter (smaller 7-ray variation for 
a given X-ray one). This might explain the apparent change 
of the flux-flux correlation trend between the beginning and 
end of the week-long campaign. 

However, it is worth noting that Krennrich et al. (2003) find 
that during the "flare" of March 25, a flux variation larger than 
a factor of 2 does not seem to be accompanied by any spec- 
tral change. Given the characteristics of the spectra, namely 
the fact that the IC peak at most moved marginally within the 
observed band, this achromaticity can not be convincingly as- 
cribed to the fact that Whipple was observing the lower energy 
shoulder of the IC peak. The possibility that it is intrinsic has 
to be contemplated. 

Therefore the possible change of the X-ray /TeV flux corre- 
lation may also be attributed to some intrinsic effect, possibly 
related to the longer term increase of luminosity. 

3.5. X-ray vs. TeV spectra and spectral energy distributions 
3.5.1. Intranight X-rayAy-ray spectra pairs 

Besides the unprecedented quality of the X-ray and 7-ray 
light curves that we have illustrated and discussed in the pre- 
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ceding sections, the March 2001 dataset affords us a unique 
opportunity of following the spectral evolution itself, with 
a time resolution that allows meaningful intra-flare analysis. 
Detailed SED-snapshot and time dependent modeling analy- 
ses are beyond the scope of this paper and will be presented in 
a forthcoming publication. Here we present the subset of X- 
ray/7-ray spectra for the March 19 event (Whipple), and for 
the March 21/22, 22/23 flares (HEGRA, presented by Aharo- 
nian et al. 2002). 

A summary "gallery" of the pairings of X-ray and 7-ray 
spectra for these flares is shown in Figure 9. For reference we 
plotted also some historical observations. 

It is worth noting that this gallery does not include the high- 
est luminosity X-ray states, nor in general (i.e. irrespective of 
simultaneous 7-ray data), neither among the intervals match- 
ing TeV observations. On the other hand, the peak of the 
March 19 flare does constitute the most luminous TeV spec- 
trum of the 2001 campaign, and in fact it matches the spec- 
trum and luminosity of the most intense flare ever recorded 
for Mrk 421, fliat of May 7, 1996 (Zweerink et al. 1997). 

For ease of comparison we prefer to adopt the same axis 
scales for X-ray and 7-ray, and this makes the variability of 
the RossiXTE spectra not as easily noticeable as that of Whip- 
ple/HEGRA spectra. Nevertheless the level of variability can 
be appreciated by comparison with the reference historical 
spectra. 

We would like to highlight a few observational findings. 
The peaks of the synchrotron and IC components never cross 
into the telescopes bandpasses, despite the relatively large lu- 
minosity variations. Increases in brightness are accompanied 
by significant spectral hardening, but there is no compelling 
sign that this is also accompanied/due by a shift of the SED 
peak energies. Among the data presented here, the only in- 
stances when the synchrotron peak might be/is directly de- 
tected are the spectra for the March 22/23 flare. It indeed 
seems that the high energy tails move between hard and soft 
states as if pivoting with respect to unobserved lower energy 
parts of the spectrum, possibly the synchrotron or IC peaks. 
This is suggested by the observation that in most cases the 
lowest energy data point in successive spectra are approxi- 
mately at the same level, whereas we would expect some "up- 
ward shift" in both the case of variations due to a change of 
energy the SED peak, and the case of an overall increase of 
luminosity around the SED peak. 

3.5.2. Spectral Energy Distributions 

In Figures 10 we show selected simultaneous X-ray and 7- 
ray spectra for the March 2001 campaign, together with a col- 
lection of historical multiwavelength data (see Figure caption 

for details). 

In particular Figure 10a shows the data for the peak of the 
March 19 outburst, and Figure 10b the pre-flare interval for 
March 22/23. These two are quite representative of a bright 
and hard, and a fainter and soft cases. We tried to model 
this sparse SEDs with a single zone homogeneous SSC model 
(e.g. Ghiselhni et al. 1998), and example fits are plotted along 
with the data. 

Although the simultaneous data coverage is limited to op- 
tical flux and the X-ray and 7-ray spectra, a coarse search 
of the parameters space for a good SSC model fit showed 
that the constraint are nonetheless very strong. This is true 
even though we made no attempt at taking into account self- 



consistently the abundant information available "along the 
time axis", such as the time resolved spectral variability. One 
general difficulty encountered while fitting the SSC model, 
is that the TeV spectra are typically harder than what can be 
predicted. As we illustrate in the section §3.5.4, this is in part 
due to the effect of the Klein-Nishina (K-N) decrease of the 
scattering efficiency, canceling the contribution from the self- 
Compton of the electrons and photons emitting/emitted above 
the synchrotron peak. 

3.5.3. B— 5 diagnostic plane 

Since we can estimate the energies and luminosities of the 
synchrotron and IC peaks with reasonable accuracy, in the 
context of a single zone SSC model we can draw the locus 
allowed by a given SED in the B - (5 parameter space (e.g. 
Tavecchio et al. 1998). Besides this primary piece of infor- 
mation, measurements or estimates of several other quantities 
(and their combinations) can be exploited to set additional 
constraint on the B -6 relationship. These include for in- 
stance peak luminosities, cooUng times, variabiUty timescales 
(or source size), intraband time lags. 

An example is shown in Figure 11, for March 22/23. The 
gray band represents the constraints set by our estimate of the 
peak positions. The peak luminosities yield a few additional 
lines in the B-5 plane, in particular the particle-magnetic field 
equipartition. A requirement on the cooling time of the peak- 
emitting particles, for instance to be shorter than the "typical" 
variability timescale (e.g. lOks), translates into an excluded 
wedge in the lower right part of the plane. The two differ- 
ent lines, meeting at the equipartition line, correspond to syn- 
chrotron or IC dominated cooling regimes. 

Detection or an upper limit on the value of intraband X-ray 
lags also sets a lower boundary to the allowed region in the di- 
agram. In Figure 1 1 we draw the limit for a hypothetical 2ks 
lag within the PCA bandpass. Shorter lags, or upper limits 
on them, move this line upwards. It is worth noting that the 
detailed cross correlation analysis of the X-ray dataset (F08) 
does not yield any reliable intraband lag detection. In partic- 
ular, no lags have been found in the analysis of all short, sin- 
gle orbit, subsets with significant variability features (a few 
dozen), and in most cases the upper limits is of the order or 
a few hundred seconds, ~ 200 -300 s (the corresponding line 
in Figure 1 1 would be about 0.5 decades, x3, higher). 

Finally, we can draw in the B-5 plane the dividing line 
between the Thomson and Klein-Nishina scattering regimes, 
for the SED peaks. 

One of the largest sources of uncertainty for the determina- 
tion of allowed region in the B-d plane is the position of the 
IC peak, because 



for scattering in Klein-Nishina regime. This also means that 
any consideration based on this diagnostic plane is subject to 
the uncertainty about the details of the TeV photon absorption 
by the diffuse infrared background. The models shown here 
include the effect of the IR absorption, following the "low 
intensity" model prescription of Stecker & De Jager (1998). 
For our limited modeling purpose the exact choice of IR back- 
ground absorption model is not critical. 

Figure 1 1 shows the locii and limits obtained for a SED sim- 
ilar to that on the right panels of Figure 10 (March 22/23), 
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Fig. 9. — Gallery of RossiXTE, Whipple/HEGRA spectra pairs, in blue symbols. Time elapses left to right, top to bottom. The top six panels 
refer to March 19 (Whipple): approximate times are (UTC) 05:39, 07:04, 07:34, 08:33 (flare top), 08:59, 09:30. The bottom two pairs are 
for March 21/22 and 22/23, with HEGRA data (Aharonian et al. 2002, preflare and flare). For reference we also plot: i) The simultaneous 
observations of the May 1994 reported by Macomb et al. (1995) (maroon, 3-point X-ray spectrum and single flux point at 0.4 TeV, at ~ 10"'"). 

ii) The highest state observed in 1996 (Zweerink et al. 1997) (orange. X-ray power law from ASCA, and 0.4-10 TeV spectrum by Whipple). 

iii) Denser-points light-gray 0.7-10 keV X-ray spectra are (bottom to top) lowest and highest state during BeppoSAX 1998 campaign (Fossati 
et al. 2000b), and the highest BeppoSAX 2000 state (Fossati et al., in preparation). 



X-RAY/TeV VARIABILITY OF Mrk 421 



17 



\og(u) [Hz] 

10 12 14 16 18 20 22 24 26 28 



\og{u) [Hz] 

10 12 14 16 18 20 22 24 26 28 




10-10 



1 10 
Energy [keV] 



1 10 
Energy [TeV] 



Energy [keV] 



1 10 
Energy [TeV] 



Fig. 10. — Spectral Energy Distributions for two epochs during the March 2001 campaign. Left, time around the peak of the March 19 
flare. Right, the "low" -preflare- state observed by HEGRA on Match 22/23 (HEGRA spectra from Aharonian et al. 2002). Simultaneous 
2001 data are shown in blue. For the optical we plot the highest and lowest fluxes observed during the campaign. The sparser, connected, 
light blue, 3-20 keV X-ray spectra show the highest and lowest X-ray observed states. For reference, we plot in gray multiwavelength data 
from a collection of historical data from NED and Macomb et al. (1995). The maroon points are the simultaneous observations for May 1994 
(Macomb et al. 1995). For the description of some of the X-ray and 7-ray data see also Figure 9. The solid red lines represent "fits" with a 
simple one-zone homogeneous SSC model with S ~ 0. 1 -0. 15 G, 5 = 20-25, Rb\ob = 10'* cm. The dashed green SED models are for "extreme" 
cases, with B~IG,S= 100, Rb\ob = 0.5 - 1 x 10'^ cm. 



for which we obtained a satisfactory SSC model fit. The 
relevant observational parameters are reported in the plot. 
All constraints are satisfied by a model having a magnetic 
field of Z? ~ 0.15G, a Doppler factor S ~ 20, a blob size of 
R ~ 10'^ cm, i.e. within the range regarded as standard in 
SSC modeling (see circle in Figure 11, solid line model in 
Figure 10b). 

This analysis shows that for this choice of parameters the 
scattering producing the TeV emission occurs in the Klein- 
Nishina regime. It is in principle possible to shift the "sweet 
spot" in the upper corner of the diagram, into the Thom- 
son regime region, by adopting a much smaller source size 
(/? ~ 5 X 10'^ cm), and in turn B ~ 2 G and 5 ~ 100 (square in 
Figure 11), and indeed a similarly satisfactory SSC fit to the 
snapshot SED can be obtained (dashed model in Figure 10b). 

A similar analysis was performed for the data of the March 
19 flare peak, shown in Figure 10a, and also in this case the 
SED could be fit both with "standard" (B = 0.1 G, (5 ~ 20, 
R = lO"" cm) "extreme" (B=l.OG,S= 100, R = 10"* cm) pa- 
rameters (coiTesponding SEDs are shown as solid and dashed 
lines in Figure 10.) 

Other considerations can help to discriminate between these 
scenarios, for instance arguments concerning time variability 
properties or the viability of having such an extreme Doppler 
factor and blob size. 

3.5.4. TeV spectral decomposition analysis 

In order to try to understand the observed correlated vari- 
ability between X-ray and TeV fluxes, it is interesting to take 
a deeper look at the composition of the emission in the TeV 



band, in terms of which electrons and seed photons contribute 
to the flux at different energies. This analysis is somewhat 
model dependent, and we are only showing it for the standard 
parameter choice introduced before. 

The idea is similar to the treatment discussed by Tavecchio 
et al. (1998), where they split the IC component in four com- 
ponents, produced by the combinations of electron and syn- 
chrotron photons below (L) and above (H) the synchrotron 
peak (for electrons the split is done at the energy mapped to 
this latter). For TeV blazars the conditions are such that the 
component H,H (electrons and photons both above the peak) 
is strongly depressed, and becomes negligible. The same 
holds true for the L,H component (Tavecchio et al. 1998). 

The same approach can be extended to an arbitrary split of 
the primary components, aiming at identifying more precisely 
the origin of the electrons and photons. 

In Figure 12 we show the IC peak for the same model on 
which we have focused above. We adopt as breaking point 
for the electron spectrum the energy corresponding to an ob- 
served synchrotron emission at 2 keV, because this is just be- 
low the RossiXTE/PCA bandpass and this splitting allows 
us to divide the electrons between those whose synchrotron 
we observe and those we don't. We consider all electrons, 
but restrict the allowed seed photons to those in the range 
5-500eV observer's frame. Different line types coiTespond 
to different ranges of seed photons within this band. The long- 
dash black lines show the full emission by 5-500 eV photons 
and electrons below and above 2 keV, with its decomposition 
in the two contributions. It is clear that most of the emis- 
sion at around 1 TeV is accounted for by these photons scat- 
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Fig. 11. — B-S plane for a set of parameters (fpcak.x, fpcak,7, ipeak.x, 
ipeak,7, representative of the 2001 campaign. The gray locus cross- 
ing the plane shows the constraint set by the synchrotron and IC 
components peak frequencies (adopted values are shown atop the 
figure box.) The long-dash steep lines show the range allowed by 
the inferred peak luminosities. The dot-dashed steep lines marks the 
equipartition between Ub and t/rad- There are two sets of long-dash 
and dot-dashed lines: the left ones correspond to a standard case of 
R = lO'* cm. Those in the right part of the figure refer to the case 
of R = 5 X lO" cm. The dotted regions are meant to represent the 
effect of a factor of 3 uncertainty on the main parameters. The thick 
black line is the approximate analytical boundary between Thom- 
son and Klein-Nishina scattering regimes. The gray dot-dashed lines 
marks the combination of parameters yielding a lOks cooling time 
for electrons emitting the synchrotron peak. The bottom-right corner 
wedge delimited by them is non-consistent with this imposed (puta- 
tive) limit. The short-dash black line is the lower bound allowed by a 
hypothetical X-ray intraband lag of 2 ks. The black circle marks the 
parameter choice for the SED model drawn solid red in Figures 10b. 
The black square instead marks a possible choice of parameters in 
the Thomson regime region of the B-S plane, dashed green SED in 
Figures 10b. 

tered by sub-keV electrons. The gray short-dash and dot-dash 
lines show a further split of the photons in the 5 - 50 eV and 
50-500eV bands, with the former contributing about 60% of 
the total at around 1 TeV. 

Because of the Klein-Nishina effect, there is no emission 
by self-Compton by electrons observed in synchrotron in the 
X-ray band. Moreover, the (non-self) IC contribution by X- 
ray-observed electrons is weak and at high enough energy that 
we can regard is as iiTelevant to the effect of the observed TeV 
rate variability. This latter is dominated by a lower energy 
section of the bandpass, around ~ 1 keV. 

The K-N depression of the scattering of the higher energy 
photons and electrons cuts off the H,H and H,L contributions, 
and makes the TeV spectrum always steeper than the X-ray 
one, at least in the more minimalist scenarios. Indeed in first 
approximation the power law of the high-energy tail of the IC 
component tends to a value q:kn — 2a2-a[, where ai and a2 
are the spectra indices below and above the synchrotron peak 
(Tavecchio et al. 1998). This is always steeper than 0:2, and 
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Fig. 12. — Modeling of the inverse Compton peak for the SSC 
model shown solid red in Figure 10b. The IC emission is split in dif- 
ferent components a la Tavecchio et al. (1998). The long-dash black 
lines show the emission obtained by restricting the synchrotron seed 
photons to the 5-500 eV range. It is split in two by considering sepa- 
rately electrons emitting in synchrotron below and above 2 keV: the 
lower energy section is thus coming from 5-500 eV photons scat- 
tered by electrons up to 2 keV. This latter is then further split into the 
contributions from 5-50 eV (light gray short-dashed line) and 50- 
500 eV photons (darker gray dot-dashed). The thick solid line is the 
full IC component, i.e. including all seed photons contributions. The 
vertical dotted lines mark the 0.4^ TeV band. 

for fiducial values of a\ and aa it is so by Aa > 0.5. 

In the more extreme scenario contrived to shift the scat- 
tering into the Thomson regime, the IC peak composition is 
indeed different. In particular, as expected, there is a more 
even contribution by the three of the four components. The 
H,H components is still negligible. Photons up to 2-3keV 
are effectively scattered, and there is a self-Compton contri- 
bution to the emission at ~ 1 keV. Therefore, a much more 
direct connection between the two observed bands is afforded 
by the more exotic scenario, and in turn the observed coiTe- 
lated variability would be more readily explained. 

4. CONCLUSIONS 

The correlation between the variations in the X-ray and TeV 
bands is confirmed with unprecedented detail, supporting the 
idea that the same electron distribution, in the same physical 
region, is responsible for the emission in both energy bands. 
However the details of these findings pose a serious challenge 
to the emission models. Here we would like to sketch a few 
selected outstanding issues raised by the correlated variabil- 
ity. We refer to a forthcoming paper for an in-depth analysis 
comprising more extensive modeling. 

• If modeled within the realm of standard values for mag- 
netic field, Doppler factor, source size, the IC scattering 
responsible for the observed TeV emission occurs in the 
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Klein-Nishina regime. This means that with X-ray and 
7-ray observations, although we seem to be observing re- 
gions of the spectrum that are very similar to each other for 
what concerns their position with respect to the SED peaks, 
we are not tracking the evolution of the same electrons (and 
photons). The extent of the phase and amplitude correla- 
tion of the X-ray and 7-ray variations is however remark- 
able, and this sets broad constraints on the characteristics 
of the processes responsible/governing the variability, e.g. 
acceleration/injection of particles, dominant cooling cause. 

In this context, the observation in K-N regime of a 
quadratic relationship between synchrotron and IC varia- 
tions (which would be naturally produced in the Thomson 
regime because of the effectiveness of self-Compton) con- 
strains the electron spectrum variations to occur over an 
energy band broad enough to affect also the IC seed pho- 
tons. Moreover, this variation must be essentially "achro- 
matic" (i.e. just a change in normalization), otherwise the 
extra energy-dependent factor would produce an observ- 
able effect. 

The observation that the flux-flux path of the better ob- 
served flares decay foUows closely the bursting path in- 
troduces a further complication. If the flare decay is gov- 
erned by the cooling of the emitting electrons we do not 
expect the quadratic relationship to hold during the decay- 
ing phase. In fact, given the energy-dependent nature of 

-1/2 

synchrotron (and IC) coolmg, with TcooI ^ E^^^' , the 50 eV 
seed photons cool on a longer timescale, e.g. ^10 times 
longer than the timescale for photons observed at > 3 keV. 
The possibility that also the electrons contributing to the 
bulk of the TeV emission have lower energy than those 
observed in X-ray, compounds the problem. This means 
that during the flare decay the X-ray and 7-ray brightnesses 
should follow something like a linear relationship, because 
the IC (TeV) emission will just reflect the evolution of the 
electron spectrum, scattering a "steady" seed photon field. 

Plain radiative cooling does not seem to match these obser- 
vations. A viable mechanisms explaining the flares evolu- 
tion should allow the concurrent cooling of a broad portion 
of the electron distribution. On the other hand, brightness 
variations are accompanied by large spectral changes, and 
in most cases they are very suggestive of acceleration -or 
injection- of the higher energy end of the electron popula- 
tion. 

Another recurring discrepancy between data and simple 
one-zone SSC modeling is that of the TeV spectral shape, 
which is often harder than model predictions. As we il- 
lustrated in the previous section the Klein-Nishina effect 
plays an important role in this respect. A more careful 
analysis is warranted, but it is worth noting that one alter- 
native option for addressing this problem is that of consid- 
ering the effect of additional IC components, off photons 
external to the blob. Blazejowski et al. (2005) showed that 
a multi-component model seems to be required to fit the 
2003-2004 observations, and it might mitigate the discrep- 
ancy in the TeV spectrum. Ghisellini et al. (2005) discuss 
the effect of including the effect of the radiation emitted 
by the putative lower Lorentz factor outer layer of the jet, 
namely as source of additional seed photons for IC. 



• Exotic scenarios could address some of these issues, 
namely by allowing the IC scattering to occur in the Thom- 
son regime, and hence the self-Compton to be effective, 
thus reinstating the close relationship between the photons 
and electrons tracked by X-ray and 7-ray observations. 
One such scenario that we discussed briefly would call for 
very high values of the Doppler beaming factor, that would 
reduce the intrinsic energies at play, higher magnetic field 
and very small size of the emission region. In fact recently 
there has been some interest for less conventional model- 
ing of TeV blazars (e.g. Krawczynski et al. 2001; Rebillot 
et al. 2006 for Mrk 421). However, high Doppler factor 
scenarios raise a series of new issues, or re-open some that 
have been settled for the more traditional model. First of 
all, the fact that the beaming cone of the radiation emit- 
ted by such a fast blob is going to be much narrower has 
to be reconciled with the population statistics of blazars 
and radio-galaxies, their putative parent population (Urry 
& Padovani 1995). One way of doing it would be to imag- 
ine that the jet comprises a very large number of small 
high Lorentz factor blobs, fanning out filling a wider cone, 
with aperture consistent with the unification statistics. This 
would constitute a quite radical change in the jet structure, 
from the current one where emission is thought to come 
from internal shocks. There are also impUcations concern- 
ing the statistical properties of the variability, which would 
likely be due to the combination of the bursting of different 
blobs, most hkely uncorrected. 

The richness and depth of the X-ray and 7-ray data of the 
March 2001 campaign presented in this paper raise the bar for 
models. The aggregate characteristics illustrated here already 
challenge the simple traditional SSC model, and the SED- 
snapshot approach. In order to answer the questions raised 
by these observations it is of paramount importance to exploit 
fully the time-axis dimension into the modeling and take a 
dynamical approach. 

The data time density and brightness (and so statistics) 
are unparalleled, enabling time resolved spectroscopy on 
timescales of the order of the physically relevant ones, hence 
allowing to model the phenomenology self-consistently min- 
imizing the need (freedom) to make assumptions as to how to 
connect spectra taken at different times. 

It is hkely that this dataset is going constitute the best 
benchmark for time dependent modeling for some time, de- 
spite the great progress made by ground based TeV atmo- 
spheric Cherenkov telescopes in the last few years, because 
of the difficulty of securing long uninterrupted observations 
with Chandra and XMM-Newton. 
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TABLE 1 

Journal of RossiXTE Observations' 



Obs. # Date(UTC) PCA HEXTE 





Start 


Stop 


Tex 

cxp 


# of GTls 


Ton 


# of GTls 


1 


18/()3/()l:21:40 


19/03/01:05:40 


16544 


6 


16544 


6 


2 


19/03/01:05:40 


19/03/01:06:11 


1472 


1 


1824 


1 


3 


19/03/01:06:55 


19/03/01:14:49 


15056 


5 


16400 


5 


4 


19/03/01:14:49 


19/03/01:16:17 


3280 


1 


3280 


1 


5 


19/03/01:18:13 


19/03/01:20:36 


5488 


2 


5984 


2 


6 


19/03/01:21:19 


20/03/01:05:19 


16464 


6 


16464 


6 


7 


20/03/01:05:19 


20/03/01:06:05 


2144 


1 


2640 


1 


8 


20/03/01:06:59 


20/03/01:08:11 


1936 


1 


2576 


1 


9 


20/03/01:08:31 


20/03/01:13:00 


8432 


3 


9360 


3 


10 


20/03/01:13:29 


20/03/01:15:26 


4560 


2 


4560 


2 


11 


20/03/01:16:20 


20/03/01:18:51 


5984 


2 


6352 


2 


12 


20/03/01:21:07 


21/03/01:05:07 


16544 


6 


16544 


6 


13 


21/03/01:05:07 


21/03/01:05:58 


2304 


1 


2832 


1 


14 


21/03/01:06:40 


21/03/01:14:32 


14720 


5 


16448 


5 


15 


21/03/01:14:32 


21/03/01:19:01 


9440 


3 


9888 


3 


16 


21/03/01:22:52 


22/03/01:02:42 


8784 


4 


8784 


3 


17 


22/03/01:03:32 


22/03/01:04:02 


1760 


1 


1760 


1 


18 


22/03/01:05:19 


22/03/01:05:48 


960 


1 


1600 


1 


19 


22/03/01:06:41 


22/03/01:12:45 


10656 


4 


12560 


4 


20 


22/03/01:13:13 


22/03/01:14:20 


2112 


1 


2112 


1 


21 


22/03/01:16:07 


22/03/01:16:46 


2280 


1 


2280 


1 


22 


22/03/01:19:19 


22/03/01:23:48 


9632 


3 


9632 


3 


23 


23/03/01:00:07 


23/03/01:03:50 


8456 


3 


8456 


3 


24 


23/03/01:06:24 


23/03/01:14:15 


14368 


5 


16368 


5 


25 


23/03/01:14:15 


23/03/01:18:25 


9512 


3 


9784 


3 


26 


23/03/01:19:06 


24/03/01:02:21 


16296 


5 


16296 


5 


27 


24/03/01:03:02 


24/03/01:03:51 


2256 


1 


2672 


1 


28 


24/03/01:04:41 


24/03/01:10:54 


9984 


4 


12784 


5 


29 


24/03/01:11:21 


24/03/01:14:54 


7744 


3 


7856 


3 


30 


24/03/01:15:45 


24/03/01:16:37 


2880 


1 


2880 


1 


31 


24/03/01:19:01 


24/03/01:21:06 


5064 


2 


5064 


2 


32 


24/03/01:22:26 


25/03/01:01:52 


7448 


3 


7448 


3 


33 


25/03/01:03:12 


25/03/01:03:42 


1152 


1 


1680 


1 


34 


25/03/01:04:55 


25/03/01:05:49 


832 


1 


1616 


2 


35 


25/03/01:06:27 


25/03/01:13:57 


13024 


5 


15216 


5 


36 


25/03/01:13:57 


25/03/01:15:33 


3296 


1 


3296 


1 


Total 


l8/03/0l:il:40 


25/03/01:15:33 


262864 


98 


2m40 


99 



^ Non deadtime corrected, and therefore strictly not an exposure time, but an on-source time. 



TABLE 2 

RossiXTE /PCA and Whipple, HEGRA Overlap Statistics 



Night # 


Date^ 


Date (MJD) 


Whipple 


Overlap Fraction 


HEGRA 


Overlap Fraction 








Exp. Time 


detailed" 


run-by-run'^ 


Exp. Time 


detailed" 


run-by-run'^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


All 


March 18-25 


51986/51993 


29h52m 


44 % 


47/ 64 (73 %) 


3244" 


54% 


58/69 (84%) 


1 


March 18/19 


51986/51987 


61104'" 


42% 


11/13 (85%) 


5" 14° 


55% 


11/11 (100%) 


2 


March 19/20 


51987/51988 


3h44m 


67 % 


7/ 8 (87%) 


4h33m 


60% 


8/10 (80%) 


3 


March 20/21 


51988/51989 


5''36'" 


45% 


9/12 (75%) 


5''30™ 


54% 


10/11 (91%) 


4 


March 21/22 


51989/51990 


4hl2'" 


44% 


7/ 9 (78%) 


4^32™ 


54% 


7/10 (70%) 


5 


March 22/23 


51990/51991 


4h4om 


30% 


5/10 (50%) 


4h37m 


57% 


8/10 (80%) 


6 


March 23/24 


51991/51992 


2^2(r 


55% 


4/ 5 (80%) 


4''35'" 


58% 


9/10 (90%) 


7 


March 24/25 


51992/51993 


3''16'" 


31 % 


4/ 7 (57%) 


3hiQm 


37% 


5/ 7 (71%) 



" HEGRA observing windows typically extend across UTC midnight. All the Whipple observing windows occur after the UTC 
midnight.'' Ratio between the actual RossiXTE /PCA on-source time during Whipple (HEGRA) runs, and the actual total Whipple 
(HEGRA) observing time.'^ Fraction of Whipple (HEGRA) «0.5 hrs runs with some overlapping RossiXTE /PCA data. The overlap 
fraction range between 7% and 100%, averaging to the niunber reported in columns (5) and (8). 
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TABLE 3 

RossiXTE /PCA and Whipple+HEGRA Overlap Statistics 



Night* 


Date'' 


Date (MJD) 


TeV 


Overlap Fraction 




('^) 




Exp. Tune 


detailed" 
(5) 


run-by-run*^ 
(6) 


All 


March 18-25 


51986/51993 


62''06"" 


49 % 


104/133 (78%) 


1 


March 18/19 


51986/51987 




48% 


22/24 (92%) 


2 


March 19/20 


51987/51988 


ghiym 


63 % 


15/18 (83 %) 


3 


March 20/21 


51988/51989 


ll''06" 


48 % 


18/23 (78 %) 


4 


March 21/22 


51989/51990 


8h44m 


49% 


14/19 (74%) 


5 


March 22/23 


51990/51991 




43 % 


13/20 (65 %) 


6 


March 23/24 


51991/51992 


6'>55"> 


57% 


13/ 15 (87 %) 


7 


March 24/25 


51992/51993 


6''26"' 


34% 


9/14 (64%) 



" HEGRA observing windows typically extend across UTC midnight.'' Ratio be- 
tween the actual RossiXTE /PCA on-source time during Whipple, HEGRA runs, 
w 0.5 hrs, and the actual total observing time."^ Fraction of Whipple, HEGRA 
runs with some overlapping RossiXTE /PCA data. The overlap fraction ranges 
between 7% and 100%, averaging to the number reported in colimm (5). 



TABLE 4 

Slope of the TeV/X-ray correlation. 



Night(s) 


#data 


PCA 2-10 keV 


PCA 2-4 keV 


PCA 4-6 keV 


PCA 6-8 keV 


PCA 9-15 keV 


HEXTE 20-60 keV 


all 


105 


0.88 ±0.07 


0.93 ±0.09 


0.89 ±0.07 


0.81 ±0.07 


0.71 ±0.06 


1.02 ±0.14 


all(dT=ld) 


7 


0.83 ±0.14 


0.89 ±0.16 


0.82 ±0.14 


0.79 ±0.1 3 


0.73 ±0.12 


1.96 ±0.35 


Individual Nights 


1=' 


11 


2.84 ±0.41 


2.29 ±1.00 


2.68 ±0.28 


2.51 ±0.21 


2.20 ±0.12 


2.47 ±0.96 


1 


22 


2.26 ±0.25 


1.98 ±0.45 


2.26 ±0.1 3 


2.12±0.16 


1.90 ±0.12 


2.73 ±0.66 


4 


14 


1.56 ±0.25 


1.31 ±0.23 


1.73 ±0.27 


1.65 ±0.25 


1.61 ±0.22 


1.44 ±0.52 


5 


14 


1.67 ±0.1 6 


1.81 ±0.17 


l,6.T±0,17 


1.61 ±0.15 


1.44±0.I4 


1,72 ±0,63 


Selected subsets of Consecutive Nights 


l='+2 


26 


2.16±0.21 


1.97±0.38 


2.12±0.17 


2.00 ±0.15 


1.82±0.12 


2.00 ±0.61 


1+2 


37 


2.04 ±0.16 


1.88 ±0.28 


2.03 ±0.14 


1.92 ±0.12 


1.73±0.10 


2.12±0.52 


4+5 


28 


1.74 ±0.15 


1.57 ±0.23 


1.78 ±0.15 


1.68 ±0.1 2 


1.45 ±0.10 


1.77 ±0.38 


6+7 


22 


1.00 ±0.14 


0.81 ±0.15 


1.03 ±0.14 


1.01 ±0.13 


0.98 ±0.12 


1.12±0.19 


5+6+7 


36 


0.94 ±0.09 


0.84 ±0.10 


0.98 ±0.09 


0.95 ±0.09 


0.91 ±0.08 


1.21 ±0.16 


Subsets by brightness level 


P + 2 


26 


2.16±0.21 


1.97 ±0.38 


2.12±0.17 


2.00 ±0.15 


1.82 ±0.12 


2.00 ±0.61 


6 + 7 


22 


1.00±0.14 


0.81 ±0.15 


1.03 ±0.14 


1.01 ±0.13 


0.98 ±0.12 


1.12±0.19 


3 + 6 + 7 


40 


0.99±0.14 


0.82±0.14 


1.04 ±0.14 


0.98±0.14 


0.90 ±0.13 


0.74 ±0.19 



^ Whipple (flare) only. 



